OPEN 



Heredity (2014) 112, 89-98 

© 2014 Macmillan Publishers Limited All rights reserved 0018-067X/14 



www.nature.com/hdy 



ORIGINAL ARTICLE 

Evolution of body shape in sympatric versus non-sympatric 
Tropheus populations of Lake Tanganyika 

M Kerschbaumer 1 , P Mitteroecker 2 and C Sturmbauer 1 

Allopatric speciation often yields ecologically equivalent sister species, so that their secondary admixis enforces competition. 
The shores of Lake Tanganyika harbor about 120 distinct populations of the cichlid genus Tropheus, but only some are 
sympatric. When alone, Tropheus occupies a relatively broad depth zone, but in sympatry, fish segregate by depth. To assess the 
effects of competition, we studied the partial co-occurrence of Tropheus moorii 'Kaiser' and 'Kirschfleck' with Tropheus polli. 
A previous study demonstrated via standardized breeding experiments that some observed differences between Tropheus 
'Kaiser' living alone and in sympatry with T polli have a genetic basis despite large-scale phenotypic plasticity. Using geometric 
morphometries and neutral genetic markers, we now investigated whether sympatric populations differ consistently in body 
shape from populations living alone and if the differences are adaptive. We found significant differences in mean shape 
between non-sympatric and sympatric populations, whereas all sympatric populations of both color morphs clustered together in 
shape space. Sympatric populations had a relatively smaller head, smaller eyes and a more anterior insertion of the pectoral fin 
than non-sympatric populations. Genetically, however, non-sympatric and sympatric 'Kaiser' populations clustered together to 
the exclusion of 'Kirschfleck'. Genetic distances, but not morphological distances, were correlated with geographic distances. 
Within- and between-population covariance matrices for T moorii populations deviated from proportionality. It is thus likely that 
natural selection acts on both phenotypic plasticity and heritable traits and that both factors contribute to the observed shape 
differences. The consistency of the pattern in five populations suggests ecological character displacement. 
Heredity (2014) 112, 89-98; doi:10.1038/hdy.2013.78; published online 25 September 2013 
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INTRODUCTION 

Phenotypic variation among individuals is produced by the combined 
effect of phenotypic plasticity and genetically based differences. 
Alternative phenotypes are subject to natural selection, no matter 
the source, leading to adaptive evolution at the population level (Via 
and Lande, 1985). Thus, variation in morphology constitutes a 
co-gradient, with environmental effects on phenotypic expression 
reinforcing or counteracting genetic differences between populations 
(Marcil et al, 2006). Novel traits can originate by environmental 
induction as well as mutation, to subsequently undergo selection and 
genetic accommodation. It has been argued that adaptive innovation 
owing to developmental plasticity may have greater evolutionary 
potential than mutationally induced ones (West-Eberhard, 2005a). In 
this respect, a timely theory of adaptive evolution must recognize the 
importance of phenotypic accommodation, that is, the refinement of 
the adaptive phenotype through subsequent mutational changes 
(Suzuki and Nijhout, 2006), in addition to selection operating on 
novel features arising from random variation owing to mutation 
(West-Eberhard, 2005b). The factors promoting or limiting plasticity 
in a given population, however, might be difficult to demonstrate 
(DeWitt et al, 1998), and it is especially challenging to assess the 
relative contribution of plasticity versus genetically based differences, 
when natural populations are compared. 

The cichlid fishes of the Great East African lakes have become a 
paradigm for explosive speciation and adaptive radiation. Two key 



innovations are thought to be responsible for the cichlid's success in 
populating the lake via adaptive radiation: their second set of jaws 
decoupled from the oral jaws (Liem, 1973) and their highly 
specialized reproductive behavior (Crapon de Caprona, 1986). These 
behavioral and morphological characteristics enable them to utilize a 
variety of specific trophic niches faster than any other fish group, so 
that complexly intertwined species communities have formed. Aside 
from geological processes, lake-level changes were argued to have 
significantly contributed to the diversification of cichlid fishes 
(Sturmbauer and Meyer, 1992; Cohen et al, 1993, 1997; Ruber 
et al, 1998; Baric et al, 2003; Verheyen et al, 2003; Sturmbauer et al, 
2005; Duftner et al, 2007), by generating repeated cycles of popula- 
tion subdivision, geographic isolation and divergence, followed by 
secondary admixis. The effect multiplied all along the lake shore 
toward a gigantic scenario of allopatric speciation, termed 'species 
pump' (Rossiter, 1995). The connection between ecomorphological 
divergence and speciation is a central issue in adaptive radiation, and 
the role of competition in driving the process has long been suggested 
(Fryer, 1959; Coulter, 1994). Competition for trophic resources has 
led to habitat partitioning and the establishment and maintenance of 
distinct trophic morphologies in the context of adaptive radiation 
(Albertson, 2008) and ecological speciation (Rundle and Nosil, 2005). 

The Lake Tanganyika cichlid fish genus Tropheus, with six nominal 
species and about 120 mostly allopatric 'color morphs' (Poll, 1986; 
Schupke, 2003), is an ideal study system to target allopatric divergence 
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and speciation at a mature stage of adaptive radiation. Despite great 
evolutionary age and a large variety of color patterns within this 
genus, morphology and size have been reported to be highly 
constrained (Sturmbauer and Meyer, 1992), owing to the fact that 
all populations and sister species fill the same trophic niche in the 
littoral species community. However, recent geometric morphometric 
studies have identified significant mean shape differences among 
several Tropheus populations, which in part might be adaptive and in 
part may be owing to neutral drift (Maderbacher et al, 2008; Postl 
et al, 2008; Herler et al, 2010). The evolutionary history and the 
phylogeography of the genus Tropheus has been investigated in a series 
of molecular genetic studies (for example, Sturmbauer and Meyer, 
1992; Sturmbauer et al, 1997; Baric et al, 2003; Sturmbauer et al, 
2005; Egger et al, 2007; Koblmiiller et al, 2008). The present 
distribution of genetic lineages shows a great degree of overlap, 
especially at the borders of the three lake basins, corroborating the 
important role of lake-level changes for dispersal and admixis 
(Sturmbauer et al, 2005). At such border regions, major lineages 
have come into secondary contact, and they have either hybridized or 
remained reproductively isolated, so that two or more species of 
Tropheus sometimes live in sympatry. 

Tropheus lives at rock and cobble shores where it inhabits a wide 
range of water depths, sometimes down to 40 m, but its highest 
density is between the depths of 0.5 and 5 m (Kohda and Yanagisawa, 
1992; Sturmbauer et al, 2008). The genus is highly specialized and 
occupies the trophic niche of mainly browsing on filamentous algae, 
whereas other genera utilize different resources in this complex littoral 
species community (Sturmbauer et al, 1992; Konings, 1998). In the 
littoral zone of freshwater lakes, several ecological factors vary strongly 
with water depth. Tropheus moorii prefers shallow water, presumably 
because of the higher algal productivity, warmer temperatures and 
lower mortality risk caused by pelagic predatory fish. Only very large 
adults tend to use slightly deeper water (>5m). 

In this paper, we have studied the effects of co-occurrence of two 
Tropheus species in the central eastern section of Lake Tanganyika: 
T. moorii (color morphs 'Kaiser' and 'Kirschfleck') and Tropheus polli. 
We investigated whether the sympatric Tropheus populations differ in 
body shape from the non-sympatric populations (that is, whether 
those populations that coexist with another Tropheus differ from 
those who do not). Thereby, we aimed to test the hypothesis that the 
coexistence of two ecologically (almost) equivalent species enforces 
competition for the available resources. Such competitive interactions 
might lead to spatial segregation and ecomorphological divergence. 
Morphological and behavioral differences among sympatric and non- 
sympatric Tropheus populations might in part be due to morpholo- 
gical plasticity and genetic variation, to be shaped by ecological or 
reproductive character displacement (Brown and Wilson, 1956; 
Pfennig et al, 2010). We studied the following natural scenario: On 
the eastern coast near the village Ikola, T. moorii 'Kaiser' lives alone 
and utilizes the entire range of preferred depth of 0 to about 5 m 
(Figure 1). Further north, T moorii 'Kaiser' occurs in sympatry 
with T polli. In this situation, T polli occupies the uppermost section 
of the rocky habitat, whereas T moorii 'Kaiser' lives in the 
deeper sections of the rocky littoral zone between depths of about 
3-5 m (Schupke, 2003; Sturmbauer, personal observations). However, 
the depth separation is not complete, and the two species show some 
overlap in their distributions. Our study also included another color 
morph, T. moorii 'Kirschfleck', which lives slightly north of T moorii 
'Kaiser' on the shores of the Mahale Mountains. Likewise, it shares the 
habitat with T polli and shows the same depth segregation as 
T moorii 'Kaiser'. T moorii 'Kirschfleck' is similar in color to 



T moorii 'Kaiser', in that T 'Kirschfleck' shows two red blotches on 
the body flanks, whereas T. 'Kaiser' shows a wide yellow band, and 
both color morphs share a relatively recent common ancestry (Egger 
et al, 2007). Note, however, that even populations assigned to the 
same color morph differ slightly from each other. Aquarium cross- 
breeding experiments confirmed assortative mating and reproductive 
isolation of T polli from T moorii 'Kaiser' and T moorii 
'Kirschfleck' (Sturmbauer, unpublished data). T. moorii 'Kaiser' and 
'Kirschfleck' readily form hybrids in captivity (Toby Veall, personal 
communication) . 

This study builds on a previous study in which we produced and 
analyzed V Y offspring of four Tropheus color morphs in a standardized 
pond environment and F : hybrids between two sets of Tropheus color 
morphs (Kerschbaumer et al, 2011). In addition to T moorii 'Mbita' 
and 'Nakaku', we bred and hybridized two populations of T moorii 
'Kaiser' — one living alone and one in sympatry with T. polli — which 
are also used in the present analysis. We found that the shift to a 
standardized pond environment not only induced marked phenotypic 
plasticity but also, at the same time, left a morphological distinctness 
among the populations intact in that the F 2 offspring of all four pond- 
raised Tropheus populations could be equally differentiated. Thus, 
despite the clear signal of phenotypic plasticity, a genetic basis for the 
observed morphometric differences between the populations exists. 
The amount of shape change due to phenotypic plasticity exceeded 
the differences among these populations by a factor of 2.4. We also 
demonstrated a heritable basis for particular morphological features 
differing among various Tropheus populations (Koch et al, 2012). 

To address the question of whether these differences have an 
adaptive background, we applied geometric morphometric methods 
in combination with an analysis of neutral genetic markers. We 
further related morphological and genetic differences to the geo- 
graphic distribution of the populations. 

MATERIALS AND METHODS 

Study populations and sampling 

During 2005-2009, about 800 individuals of two color morphs and two species 
of the genus Tropheus were collected from eight locations on the eastern coast 
of Lake Tanganyika (Figure 1). At three locations, one T. moorii population 
occupied the full range of the preferred habitat (termed 'non-sympatric 
Tropheus 'Kaiser"), and at five locations, one population shared its habitat with 
a sister species, T. polli (termed 'sympatric Tropheus 'Kaiser" and 'sympatric 
Tropheus 'Kirschfleck"). T. moorii 'Kaiser', which is also named Tropheus 'Ikola' 
in the aquarium trade, has an overall black body with a bold yellow band 
across the center of the body. The basic coloration of T. moorii 'Kirschfleck' is 
black, with two prominent red blotches on its body flanks. Both color morphs 
show a cline in hue along their distribution range. T. polli is characterized by a 
striped pattern (females and juveniles) or uniform bluish-grey color (mature 
males) and a unique deeply forked caudal fin. Sample sizes, names and 
geographic coordinates of the sampling sites are listed in Table 1. For this 
study, only adult fish were considered. Digital images of anesthetized speci- 
mens were obtained using a common flatbed scanner (Herler et al, 2007). 
A small fin clip was taken for genetic analysis. 

Morphometric analysis 

External body form was quantified by digitizing the Cartesian coordinates of 
18 landmarks (Figure 2) using TpsDig 2.10 (Rohlf, 2006). The landmark 
configurations were superimposed using a generalized Procrustes analysis 
(Rohlf and Slice, 1990; Mitteroecker and Gunz, 2009) and projected into 
Procrustes tangent space (Rohlf, 1999). The empirical distribution of popula- 
tion mean shapes was assessed using a principal component (PC) analysis of 
the resulting shape coordinates. PCs and group mean differences were 
visualized using thin-plate spline (TPS) deformation grids (Bookstein, 1991). 
Total within-population shape variance was computed as the trace of the 
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Figure 1 Sampling localities at Lake Tanganyika, (a) Tropheus moorii 'Kirschfleck' in sympatry with Tropheus polli near Mahale (KFS1) and Mabilibili 
(KFS2); (b) T. moorii 'Kaiser' in sympatry with T. polli from the south of Isonga (IKS3), north of Kekese (IKS4) and at Kekese (IKS5); sympatric T. polli 
(TPS1-TPS5); (c) Non-sympatric T. moorii 'Kaiser' living alone without a second Tropheus north of I kola (IKA1), at I kola (IKA2) and south of I kola (IKA3). 



Table 1 Summary of the sampled Tropheus species and populations 



Sampling site 


Year of 
sampling 


Code 


Living 


Coordinates 


Species 


Color morph 


Sample size 
(males/females) 


Kaiser 1 


2009 


IKA1 


Non-sympatric 


6°40'29"S, 30°20'58"E 


T. moorii 


'Kaiser' 


69 (23/46) 


Kaiser 2 


2005 


IKA2 


Non-sympatric 


6°41'27"S, 30°21'41"E 


T. moorii 


'Kaiser' 


73 (25/48) 


Kaiser 3 


2007 


IKA3 


Non-sympatric 


6°41'30"S 30°21'47"E 


T. moorii 


'Kaiser' 


141 (50/91) 


Mahale 


2009 


KFS1 


Sympatric 


6°26'46"S, 29°54'15"E 


T. moorii 


'Kirschfleck' 


54 (36/18) 




2009 


TPS1 






T. polli 




30 (13/17) 


Mabilibili 


2009 


KFS2 


Sympatric 


6°27'03"S, 29°54'51"E 


T. moorii 


'Kirschfleck' 


61 (29/32) 




2009 


TPS2 






T. polli 




28 (13/15) 


South of Isonga 


2009 


IKS3 


Sympatric 


6°30'40"S, 30°11'26"E 


T. moorii 


'Kaiser' 


66 (31/35) 




2009 


TPS3 






T. polli 




27 (15/12) 


North of Kekese 


2009 


IKS4 


Sympatric 


6°36'55"S, 30°17'35"E 


T. moorii 


'Kaiser' 


71(28/33) 




2009 


TPS4 






T. polli 




28 (12/16) 


Kekese 


2005 


IKS5 


Sympatric 


6°36'57"S, 30°17'40"E 


T. moorii 


'Kaiser' 


67 (23/44) 




2005 


TPS5 






T. polli 




48 (23/25) 



corresponding covariance matrix of shape coordinates. Population differences 
in mean shape and total variance were tested for statistical significance using 
Monte-Carlo permutation tests (Good, 2000). 

According to Lande (1979), the additive genetic between-population 
covariance matrix is expected to be proportional to the ancestral within- 
population covariance matrix under pure genetic drift. Deviations from 
proportionality may indicate evolutionary scenarios involving disruptive or 
stabilizing selection (Chapuis et al> 2008). As surrogates for the genetic 
covariance matrices, we estimated the pooled phenotypic within-population 
covariance matrix and the phenotypic between-population covariance matrix 
of the X moorii populations (Cheverud, 1988; Roff, 1997). We performed a 
likelihood ratio test of proportionality of the two covariance matrices (Mardia 
et al, 1979; Chapuis et al, 2008). In order to compute this test, the data were 
reduced to the first 10 PCs so that the covariance matrices were invertible. For 
our data, the results did not depend on the number of selected PCs. 

This quantitative genetic approach relies on several highly idealized 
assumptions (see, for example, Lande, 1979; Cheverud, 1988; Roff, 1997; 



Marroig and Cheverud, 2004; Pigliucci, 2006) and is applied to a limited 
number of populations. Hence, in addition to the likelihood ratio test, we 
explored how much the within- and between-population covariance matrices 
actually differed relative to the differences between the population covariance 
matrices. For this purpose, we performed an ordination analysis of the eight 
population covariance matrices, the pooled within-population covariance 
matrix and the between-population covariance matrix. As a distance function, 
we used the two-norm of the log relative eigenvalues, the natural metric on the 
space of covariance matrices (Mitteroecker and Bookstein, 2009). Again, this 
analysis must be based on a small set of PCs so that all matrices are invertible. 
Analyses of the first three to seven PCs yielded roughly the same ordination; 
we thus report the four PC version here. 

Furthermore, we assessed the pattern of differences among the within- and 
the between-population covariance matrices using scatter plots of the 
corresponding shape coordinates, together with landmark-wise equal 
frequency ellipses. This 'graphical decomposition' of variation allows for a 
visual inspection and localization of differences in the patterns of variance and 
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covariance (see also Herler et al, 2010). All statistical and morphometric 
analyses were performed in Mathematica 8 using routines programmed by 
Philipp Mitteroecker and Philipp Gunz. 

Genetic analysis 

Genomic DNA was extracted from ethanol-preserved fin clips via enzymatic 
digestion using proteinase K, followed by ammonium acetate and isopropanol 
precipitation (Sambrook et al, 1989). All samples were examined for genetic 
variation in six microsatellite markers: UNH130 (Lee and Kocher, 1996), 
Pzep2, Pzep3 (van Oppen et al, 1997), UME003, UME002 (Parker and 
Kornfield, 1996) and TmoM27 (Zardoya et al, 1996). PCR amplification was 
carried out in a total volume of 20 ul. PCR cocktail contained 50 ng of 
extracted total DNA, 0.025-0.2 um ul -1 of primer (forward primers were 
fluorescently labeled with HEX, 6-FAM or NED), 0.05 nmol dNTP mix 
(10 mm), 0.5 |il of BSA (lOOxBSA; New England Biolabs, Ipswich, MA, 
USA), 30 nmol of MgCl buffer (15 mm) and 0.2 ul of Taq polymerase 
(5Uul _1 ; BioTherm, GenXpress, Vienna, Austria). All PCR reactions were 
performed under the following conditions: 94 °C for 3min, followed by 30 
cycles at 92 °C for 30 s; 5 1-54 °C, 1 min; 72 °C, 1 min, followed by 72 °C for 
lOmin. PCR products were loaded on an Applied Biosystems 3130x1 sequencer 




Figure 2 Positions of the 18 digitized landmarks: (1) Anterior tip of the 
snout; (2) and (3) anterior and posterior insertion of the dorsal fin; (4) and 
(6) upper and lower insertion of caudal fin; (5) midpoint of the origin of the 
caudal fin; (7) and (8) posterior and anterior insertion of the anal fin; (9) 
insertion of the ventral fin; (10) ventral tip of cleithrum (11); most ventral 
point of the border between interoperculum and sub-operculum; (12) the 
point where preoperculum, inter-operculum, and suboperculum get in contact; 
(13) upper insertion of the pelvic fin; (14) dorsal origin of the operculum; 
(15) dorsal end of the preopercular groove; (16) and(17) most anterior and 
most posterior points of the orbit; and (18) most posterior point of the lips. 



and scored by eye using the Genemapper v.3.7 software (Applied Biosystems, 
Vienna, Austria) against an internal ABI ROX 500 size standard. We used the 
software Micro- Checker (van Oosterhout et al, 2004) to assess the potential 
presence of null alleles. The inferred genotypic information was evaluated for 
deviations from the Hardy- Weinberg equilibrium (FIS per population), and 
gene diversity, number of alleles and allelic richness were calculated using the 
software FSTAT v. 2.9.3.2 (Goudet 1995). Using the software GENEPOP 4.0 
(Raymond and Rousset, 1995; Rousset, 2010) tests for linkage disequilibrium 
were performed using the default Markov chain parameters. Among-popula- 
tion differentiation was quantified using Arlequin 3.1.1 (Excoffier et al, 2006) 
by calculating pair- wise F S t (Wright, 1951). To illustrate genetic divergence 
among populations, we performed a principal coordinate analysis based on the 
Fst distance matrix. In order to visualize the relative similarity among different 
groups of populations, a factorial correspondence analysis was performed on 
microsatellite data using Genetix v. 4.05 (Belkhir et al, 2004; data not shown). 

RESULTS 

Morphological differentiation 

Figure 3a shows a scatter plot of the first two PCs of the 13 population 
means. These two components accounted for 78% of the total shape 
variation among the population means. The sympatric 'Kaiser' 
populations were very similar to the sympatric 'Kirschfleck' popula- 
tions, whereas both differed from the non-sympatric Tropheus popula- 
tions. The mean shapes of the five T polli populations clearly differed 
from those of all T moorii 'Kaiser' and 'Kirschfleck' populations. The 
shape differences corresponding to the two PCs are visualized in 
Figures 3b and c. Despite apparent mean differences, there was an 
overlap of individual variations between the populations, both in a PC 
analysis and in a canonical variate analysis. Only T polli and T moorii 
could be separated in a two-group discriminant analysis (not shown). 

In full Procrustes shape space, the three overall mean shapes of 
sympatric 'Kaiser', non-sympatric 'Kaiser' and sympatric 'Kirschfleck' 
all differed significantly from each other (P< 0.001 for the three tests). 
Figure 4 shows TPS deformation grids between these group mean 
shapes. The non-sympatric Tropheus populations had a relatively 
enlarged head and a more posteriorly positioned pectoral fin 
compared with the sympatric populations. The differences between 
non-sympatric 'Kaiser' and sympatric 'Kirschfleck' were virtually the 
same as those depicted in Figure 2a and thus are not shown. Shape 
differences between T moorii and T polli were detected for overall 
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Figure 3 Principal component analysis of the 13 population mean shapes, (a) Scatter plot of the first two principal components (PCs), accounting for 78% 
of total shape variation among the population means. The shape differences corresponding to PC 1 and PC 2 are visualized in (b) and (c) as deformation 
grids from the grand mean shape to shapes corresponding to scores of -0.1 and 0.1 along the PCs, respectively. IKS, sympatric T. moorii 'Kaiser'; IKA, 
non-sympatric T. moorii 'Kaiser'; KFS, sympatric T. moorii 'Kirschfleck'; TPS, T. polli. 
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C T. polli. — ► T. moorii allop. T. moorii allop. — ► T. polli. 




Figure 4 Thin-plate spline visualization of mean shape differences between (a) sympatric and non-sympatric Tropheus moorii populations, (b) Tropheus polli 
and sympatric T moorii populations and (c) T polli and non-sympatric T moorii populations. All differences are linearly extrapolated by a factor of 6. 



Table 2 Total within-population shape variance (scaled by 10 3 ) 
pooled over both sexes, as well as, separately, for females and for 
males 





Pooled 


Females 


Males 


T. moorii 'Kaiser' non-sympatric 


0.462 


0.457 


0.476 


T. moorii 'Kaiser' sympatric 


0.372 


0.400 


0.341 


T. moorii 'Kirschfleck' sympatric 


0.387 


0.466 


0.338 


T. polli sympatric 


0.280 


0.291 


0.271 



body proportions, pectoral fin position, as well as the relative size and 
shape of the head. 

Table 2 gives the total variance of Procrustes shape coordinates for 
the three groups of T. moorii populations and for the T. polli 
populations. Pooled over the corresponding populations and over 
both sexes, non-sympatric 'Kaiser' had a larger total within-popula- 
tion variance than both sympatric T. moorii 'Kaiser' (P< 0.001) and 
sympatric T. moorii 'Kirschfleck' (P< 0.001). The within-population 
variance of T. polli was lower than that of all T. moorii groups 
(P< 0.001). When computed separately for the two sexes, this pattern 
was more pronounced in males than in females. 

The maximum likelihood test indicates that the between-popula- 
tion covariance matrix and the pooled within-population covariance 
matrix of all T. moorii specimens deviated significantly from 
proportionality (P< 0.001). The ordination analysis in Figure 5 shows 
that they differed quite substantially relative to the differences among 
the covariance matrices of the eight populations. The between- and 
within-population covariance structures are visualized in Figure 6 as 
scatter plots of the Procrustes shape coordinates along with equal 
frequency ellipses for each landmark. The landmark distribution 
within the populations (Figure 6a) deviated from that between the 
populations (Figure 6b) mainly in the anterior part of the head and at 
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Figure 5 Principal coordinate ordination (PCO) of the eight Tropheus moorii 
population covariance matrices, together with the pooled within-population 
covariance matrix (W) and the scaled between-population covariance matrix 
(B). Each point in this plot represents one covariance matrix, and the 
distance between the points approximates the metric presented by 
Mitteroecker and Bookstein (2009). The matrix B was scaled to be as 
proportional as possible to W by a maximum likelihood scaling factor 
(Mardia et al., 1979). 

the insertion of the pectoral fin. Note that these plots only show the 
variances of the shape coordinates along with the covariances between 
x and y coordinates of the same landmark but do not represent 
covariances between shape coordinates of different landmarks. 

Genetic differentiation 

In the two sympatric T. moorii 'Kirschfleck' populations and some 
T. polli populations, several loci showed evidence of null alleles owing 
to homozygote excess. One locus, UME003, showed evidence of null 
alleles in all but one population, so we exclude this marker from 
further analyses. A summary of microsatellite statistics is given in 
Table 3. Although genetic differentiation among non-sympatric and 
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DISCUSSION 

The relative importance of natural selection and random genetic drift 
for diversification and speciation has been the central interest of a wide 
range of empirical and theoretical studies (for example, Otte and 
Endler, 1989; Barton, 1996; Orr, 1998; Dieckmann and Doebeli, 1999; 
Schluter, 2000; Kingsolver et al, 2001; Merila and Crnokra, 2001; 
Coyne and Orr, 2004; Gavrilets, 2004; Leinonen et al, 2006). In this 
paper, we attempted to grasp the effects of natural selection on pairs of 
sympatric species by comparing the degree of phenotypic and genetic 
variation within and among selected non- sympatric Tropheus popula- 
tions with the degree of variation within sympatric populations, that 
is, those who live with another Tropheus. T polli shares its habitat with 
various Tropheus morphs along a stretch of about 100 km of the Lake 
Tanganyika shoreline, feeds on the same diet and is thus likely to be a 
strong niche competitor (Axelrod, 1977; Poll, 1986). Our population 
sample included populations of two of these color morphs, T 'Kaiser' 
and T. 'Kirschfleck'. When alone, Tropheus populations occur in a 
broad depth range, with the highest density occurring within the depth 
range of 0.5-5 m, but when they co-occur with another Tropheus, one 
always occupies the uppermost and seemingly preferred depth range, 
whereas the second (or third) Tropheus occupies deeper water, albeit 
with a zone of overlap. Sympatric populations thus inhabit a different 
and in fact more restricted habitat than the non-sympatric popula- 
tions. In our case, T polli always remains in the shallow depth zone, 
whereas the second (or third) Tropheus species moves toward deeper 
water. The movement toward deeper water is evident from the T. 
moorii 'Kaiser' populations that live alone, because these inhabit the 
entire depth range of rock habitat, typical for all non-sympatric 
Tropheus all around the lake. The same is true for Tropheus duboisi in 
the northern and northeastern part of Lake Tanganyika; T duboisi 
consistently lives in deeper water than the second (or third) Tropheus 



Table 3 Statistics for the five microsatellites of the investigated Tropheus populations 





Population 


IKA1 


IKA2 


IKA3 


IKS3 


IKS4 


IKS5 


KFS1 


KFS2 


TPS1 


TPS2 


TPS3 


TPS4 


TPS5 


All 




N 


70 


82 


141 


69 


71 


67 


56 


67 


24 


27 


30 


28 


75 


732 


UNH 130 


A 


17 


20 


18 


23 


20 


23 


16 


19 


12 


9 


8 


15 


19 


35 




AR 


11.369 


11.848 


10.552 


17.646 


14.749 


15.597 


13.452 


15.682 


12.000 


8.625 


7.524 


14.209 


13.815 


18.228 




HE 


0.785 


0.798 


0.748 


0.930 


0.889 


0.909 


0.906 


0.918 


0.845 


0.532 


0.705 


0.874 


0.819 






HO 


0.814 


0.793 


0.759 


0.899 


0.873 


0.896 


0.750 


0.582 


0.458 


0.444 


0.700 


0.607 


0.467 






Fis 


-0.030 


0.012 


-0.011 


0.041 


0.024 


0.022 


0.180 


0.372 


0.474 


0.183 


0.024 


0.322 


0.436 




Pzep3 


A 


5 


3 


5 


6 


4 


5 


4 


6 


3 


7 


8 


5 


5 


10 




AR 


4.951 


2.999 


4.724 


4.811 


3.949 


4.097 


3.421 


4.059 


3.000 


6.200 


7.150 


4.624 


3.859 


5.216 




HE 


0.727 


0.535 


0.694 


0.565 


0.629 


0.638 


0.445 


0.624 


0.618 


0.591 


0.353 


0.222 


0.621 






HO 


0.700 


0.537 


0.872 


0.681 


0.747 


0.687 


0.544 


0.868 


0.862 


0.833 


0.400 


0.207 


0.960 






Fis 


0.044 


0.004 


-0.254 


-0.199 


-0.180 


-0.069 


-0.213 


-0.383 


-0.381 


-0.397 


-0.117 


0.087 


-0.541 




Pzep2 


A 


15 


17 


17 


18 


19 


19 


15 


19 


19 


19 


17 


17 


21 


32 




AR 


11.467 


12.683 


11.605 


14.029 


14.170 


15.847 


12.955 


16.066 


18.115 


17.612 


16.067 


16.193 


17.590 


18.35 




HE 


0.857 


0.871 


0.869 


0.903 


0.896 


0.918 


0.908 


0.924 


0.929 


0.923 


0.914 


0.915 


0.933 






HO 


0.985 


0.963 


0.929 


0.971 


0.986 


0.925 


0.932 


0.955 


0.931 


0.967 


1.000 


0.966 


0.933 






Fis 


-0.142 


-0.101 


-0.065 


-0.068 


-0.093 


0.000 


-0.019 


-0.026 


0.015 


-0.031 


-0.077 


-0.038 


0.007 




TmoM27 


A 


6 


6 


6 


7 


6 


8 


6 


9 


4 


3 


3 


4 


5 


12 




AR 


4.795 


5.170 


4.653 


5.799 


5.375 


6.466 


4.989 


6.416 


3.846 


3.000 


2.993 


3.842 


3.453 


6.205 




HE 


0.615 


0.635 


0.574 


0.735 


0.755 


0.773 


0.472 


0.657 


0.501 


0.292 


0.213 


0.146 


0.151 






HO 


0.643 


0.610 


0.582 


0.826 


0.789 


0.750 


0.356 


0.338 


0.077 


0.259 


0.233 


0.154 


0.147 






Fis 


-0.039 


0.046 


-0.009 


-0.116 


-0.037 


0.038 


0.254 


0.490 


0.852 


0.129 


-0.080 


-0.036 


0.038 




UME002 


A 


3 


4 


4 


4 


6 


5 


3 


8 


3 


10 


7 


6 


8 


17 




AR 


2.570 


2.585 


2.742 


3.790 


4.974 


3.816 


2.873 


4.985 


2.800 


9.424 


6.746 


5.651 


6.208 


5.818 




HE 


0.389 


0.426 


0.425 


0.306 


0.384 


0.273 


0.201 


0.231 


0.099 


0.631 


0.385 


0.533 


0.671 






HO 


0.371 


0.366 


0.355 


0.275 


0.352 


0.250 


0.220 


0.221 


0.103 


0.379 


0.222 


0.379 


0.720 





Abbreviations: A, number of alleles; AR, allelic richness; Fis, Fis-values; HE, expected heterozygosity; HO, observed heterozygosity; N, sample size. 
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sympatric 'Kaiser' populations was significant (F$t between 0.033 and 
0.085), it did not exceed the differentiation within non-sympatric and 
sympatric populations (Table 3). The principal coordinate analysis of 
the F$t values showed a close genetic relatedness among all six 'Kaiser' 
populations, whereas the two color morphs, 'Kaiser' and 'Kirschfleck', 
were clearly distinct (Figure 6). T polli differed from both 'Kaiser' and 
'Kirschfleck', except for one population (TPS3) that clustered with the 
'Kirschfleck' populations. 

Geographic differentiation 

We correlated F$t values and Procrustes shape distances among the 
eight T moorii 'Kaiser' and T moorii 'Kirschfleck' populations 
with geographic distances among the corresponding sites (Figure 7 
and Table 4). F$t was strongly associated with geographic distance 
(r = 0.93, P< 0.001), whereas Procrustes distance was unrelated to 
geographic distance (r=0.05, P=0.79) and F ST (r= —0.15, 
P = 0A4). The same pattern became evident when analyzing the six 
T moorii 'Kaiser' populations alone (not shown). 




Figure 6 Scatter plots representing the variation of Procrustes shape 
coordinates (a) within Tropheus moorii populations and (b) between 
T. moorii population means. Equal frequency ellipses are separately drawn 
for each landmark. Under a scenario of pure evolutionary drift, the patters 
in (a) and (b) should be proportional, so that deviations from proportionality 
indicate selective forces. 
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PC0 3 




Figure 7 Principal coordinate ordination (PCO) of the F S j distances between all sampled populations. Distances between the points in this plot approximate 
the F S t distances between the populations. 



Table 4 F S t values (above diagonal) and Procrustes distances (below diagonal) between all 13 populations 





IKAl 


IKA2 


IKA3 


IKS3 


IKS4 


IKS5 


KFS1 


KFS2 


TPS1 


TPS2 


TPS3 


TPS4 


TPS5 


IKAl 




0.0410 


0.0182 


0.0776 


0.0522 


0.0656 


0.2546 


0.2286 


0.1831 


0.1885 


0.2762 


0.2126 


0.1689 


IKA2 


0.0178 




0.0155 


0.0680 


0.0329 


0.0421 


0.2217 


0.2116 


0.1787 


0.1787 


0.2376 


0.1720 


0.1770 


IKA3 


0.0159 


0.0141 




0.0848 


0.0420 


0.0560 


0.2386 


0.2203 


0.1769 


0.1842 


0.2568 


0.1902 


0.1704 


IKS3 


0.0105 


0.0240 


0.0216 




0.0172 


0.0175 


0.2293 


0.2144 


0.2134 


0.2130 


0.2465 


0.2149 


0.2062 


IKS4 


0.0127 


0.0246 


0.0246 


0.0106 




0.0009 


0.2226 


0.2087 


0.2011 


0.2012 


0.2304 


0.2091 


0.1997 


IKS5 


0.0121 


0.0143 


0.0116 


0.0175 


0.0194 




0.2307 


0.2139 


0.2147 


0.2160 


0.2398 


0.2272 


0.2158 


KFS1 


0.0094 


0.0212 


0.0187 


0.0098 


0.0158 


0.0148 




0.0466 


0.1095 


0.2213 


0.0839 


0.2098 


0.2304 


KFS2 


0.0095 


0.0237 


0.0214 


0.0075 


0.0095 


0.0154 


0.0082 




0.0896 


0.1869 


0.1295 


0.1894 


0.1968 


TPS1 


0.0154 


0.0253 


0.0242 


0.0148 


0.0167 


0.0194 


0.0163 


0.0137 




0.0716 


0.1429 


0.0968 


0.0881 


TPS2 


0.0189 


0.0296 


0.0305 


0.0184 


0.0161 


0.0260 


0.0210 


0.0173 


0.0112 




0.2021 


0.0683 


0.0930 


TPS3 


0.0183 


0.0265 


0.0289 


0.0188 


0.0170 


0.0242 


0.0203 


0.0175 


0.0113 


0.0065 




0.2455 


0.2761 


TPS4 


0.0186 


0.0292 


0.0307 


0.0182 


0.0173 


0.0255 


0.0206 


0.0173 


0.0109 


0.0064 


0.0070 




0.0553 


TPS5 


0.0205 


0.0218 


0.0272 


0.0218 


0.0214 


0.0218 


0.0214 


0.0204 


0.0152 


0.0175 


0.0132 


0.0152 





species (Brichard, 1978), as does one of the two sympatric Tropheus 
species near Namanzi and Mtosi on the southeastern shore of the lake 
(personal observations). The water depth is likely to enforce 
differential selection regimes, as the more shallow-living species has 
to cope with more bird predators and wave action than does the 
deep-living species but profits from higher algal growth, whereas the 
deep-living species has to deal more with pelagic predators, lower 
algal productivity and light transmission. One can thus argue that the 
observed depth segregation leads to differential adaptation. 

In our analysis, we found significant differences in mean shape 
between non-sympatric and sympatric Tropheus populations. More- 
over, the sympatric populations of both T moorii 'Kaiser' and 
T moorii 'Kirschfleck' clustered together in the PC analysis and only 



slightly overlapped in shape space with the non-sympatric popula- 
tions (Figure 3). Yet, an ordination of genetic similarities (F S t) 
showed a different picture: both non-sympatric and sympatric 
T moorii 'Kaiser' populations closely clustered together, to the 
exclusion of T moorii 'Kirschfleck' (Figure 6). This observation is 
congruent with a phylogeographic analysis of the genus Tropheus 
based on AFLP (amplified fragment length polymorphism) markers 
(Egger et al, 2007), in which a closer relationship between T moorii 
'Kaiser' and 'Kirschfleck' was found, contrasting earlier mtDNA-based 
results assigning them to different mtDNA lineages (Sturmbauer and 
Meyer, 1992; Sturmbauer et al, 1997; 2005). Although the genetic 
pattern may have originated to some extent from a combination of 
founder effects and drift, the incongruence between neutral genetic 
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similarities and morphological similarities is a strong indication of the 
action of differential selective processes in non-sympatric versus 
sympatric Tropheus populations. Despite the increase of genetic 
isolation among T. moorii 'Kaiser' and T. moorii 'Kirschfleck' 
populations with geographic distance, sympatric populations of both 
color morphs were morphologically alike and differed from non- 
sympatric populations. Sympatric populations had a relatively smaller 
head, smaller eyes and a more anterior insertion of the pectoral fin 
than did non-sympatric populations (Figure 4a). These features 
appear to be an adaptation to environmental features at greater 
water depth and reduced light transmission. For example, Jordan 
et al. (2008) observed a decrease of eye diameter with a decrease of 
light transmission in rock-dwelling cichlids from Lake Malawi. Also, 
Tropheus living in shallower water may need larger eyes to better 
detect bird predators. Differences in head size and eye size could also 
partly be an adaptation to (slightly) different food resources in 
shallow versus deeper water depth. 

We further found that genetic distances (Fst) between Tropheus 
populations were strongly associated with geographic distances, as 
expected, for neutral markers (Figure 8a). Morphological distances 
(Procrustes distances), by contrast, were unrelated to geographic 
distances (Figure 8b). The independence of genetic distances and 
morphological distances (Figure 8c) confirms the incongruence 
between the two ordination analyses in Figures 3 and 6. 

T polli differed both genetically and morphologically from the 
other species. However, one population (TPS3) shared some alleles 
with T moorii 'Kirschfleck' and hence differed from the other T polli 
populations in the ordination in Figure 6. The observed similarity to 
its sympatric ally is indicative of a past introgression event at this 
location. In fact, the presence of mtDNA haplotypes from two major 
lineages in T moorii 'Kirschfleck' suggest a short period of gene flow, 
possibly in the form of a reinforcement scenario upon initial 
secondary contact (Sturmbauer et al, 1997). 

Additional evidence of selective processes comes from the compar- 
ison between within-population and between-population covariance 
matrices for the studied T moorii populations. We found that the two 
phenotypic covariance matrices clearly deviated from proportionality, 
indicating evolutionary scenarios involving disruptive or stabilizing 
selection. Variation between populations (Figure 5b) differed from 
that of within populations (Figure 5a), most clearly in the landmarks 
of the head, particularly of the mouth, and of the anterior insertion of 
the pectoral fin. In these anatomical regions, non-sympatric and 
sympatric populations differed the most (Figure 4a). 

We also found that the total within-population variance (summed 
over all landmarks) was significantly smaller in the sympatric 
populations than in the non-sympatric populations. This can be a 
result of increased stabilizing selection pressure owing to food 



competition of T moorii with T polli and other fishes in the deeper 
water. Recent work on Lake Malawi rock- dwelling cichlids on 
community structure and phenotypic divergence have suggested that 
sedimentation seems to be another factor affecting foraging behavior 
structuring species and eco-types along this habitat gradient 
(Albertson, 2008; Parnell and Streelman, 2011), so niche segregation 
might as well happen among the two Tropheus. 

Two previous studies, in which offspring of sympatric and non- 
sympatric Tropheus 'Kaiser' were bred in a standardized pond environ- 
ment and compared with wild fish, revealed striking phenotypic 
plasticity induced by the sudden shift of several environmental 
parameters in the ponds. Lack of competitors and predators and 
completely calm water might be the most striking changes. However, 
the same experiment demonstrated a clear genetic basis of some of the 
observed inter-population differences in the Fj offspring (Kerschbaumer 
et al, 2011; Koch et al, 2012). One might argue that non-sympatric and 
sympatric populations live in slightly different habitats, so the observed 
mean shape differences could in part be due to the different environ- 
ments shaping ontogenetic development. Likewise, the reduced shape 
variation in sympatric populations might result from their decreased 
habitat heterogeneity, given that each species inhabits one part of the 
original depth range. As argued by West-Eberhard (2005a, b), such 
environmentally induced phenotypic differences, within the impressive 
scope of phenotypic plasticity, are likely to be the initial basis on which 
natural selection can act. The observed genetic basis of population 
differences demonstrates a contribution of random mutations, provid- 
ing the basis for genetic accommodation. Thus, phenotypic plasticity 
can promote divergence within populations, divergence via resource 
partitioning and character displacement upon secondary contact 
between species, to ultimately drive speciation and adaptive radiation 
(Pfennig et al, 2010). Such scenarios were postulated for the repeated 
formation of benthic and limnetic sticklebacks (Schluter and McPhail, 
1992). It was also argued that developmental plasticity mirrors the 
morphological differences found in three genera of spadefoot toads 
across two continents, as an indication of how ancestral plasticity may 
lead to morphological diversity driven by adaptation to broad environ- 
mental changes (Gomez-Mestre and Buchholz, 2006). 

It is thus evident that natural selection acts on both the scope of 
phenotypic plasticity and heritable traits and that both factors 
contribute to the observed shape differences among non-sympatric 
and sympatric Tropheus populations. Phenotypic plasticity is a highly 
important adaptive feature allowing individuals to react to more short- 
term environmental changes. Once natural selection begins to act 
consistently in one direction, for example, on account of behavioral 
changes after secondary admixis, genetically based changes are likely to 
emerge on top of plasticity, and this is exactly what our data suggest. 
We note that our study design does not allow us to disentangle the 
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Figure 8 Scatter plots of (a) geographic distances versus F S j values, (b) geographic distances versus Procrustes distances and (c) F S j values versus 
Procrustes distances between the eight Tropheus moorii populations. F S r is strongly associated with geographic distance (r=0.93), whereas Procrustes 
distance is unrelated to both geographic distance and F$t (r=0.05 and r— -0.15, respectively). 



Heredity 



Evolution of body shape in Tropheus 

M Kerschbaumer et al 



relative contributions of phenotypic plasticity and genetically based 
differences. The relative contributions of plasticity and genetically based 
differences can be disentangled by experimental studies where sympa- 
tric and non-sympatric populations are raised under various environ- 
mental conditions and, additionally, in a common garden design. 

It seems highly likely that the body shape differences between non- 
sympatric and sympatric T moorii populations are caused by 
ecological character displacement. Experimental data showing the 
fitness effects of competition and measures of selection on the traits of 
interest would be useful in making a stronger case. One alternative, 
reproductive character displacement, that is, the increase of isolation 
between taxa that are already good species (as opposed to reinforce- 
ment), seems less likely, given that all distinctive traits have a plausible 
ecological background. According to the six check-criteria for 
character displacement suggested by Grant (1972) and Arthur 
(1982) and considering the remarks of Schluter and McPhail 
(1992), the following evaluation can be made: 

1. Chance can be ruled out as an explanation of the pattern, given 
that other cases of sympatry also constitute the same depth 
segregation of the two entities. 

2. According to pond-breeding experiments (Kerschbaumer et al, 2011; 
Koch et al, 2012), the phenotypic differences among populations in 
sympatry and allopatry have a genetic basis, on top of plasticity. 

3. That enhanced differences between sympatric species are the 
outcome of evolutionary shifts, not simply of the inability of 
similarly sized species to coexist (for example, because of compe- 
titive exclusion or introgression), seems likely owing to the 
adaptive background of the changes, but must be investigated 
further by fitness- related studies on traits under divergent selection. 

4. As sympatric populations of T moorii consistently occupy the 
deeper section of the original depth range, and as they share 
changes in eye diameter, head proportions and the insertion of the 
pectoral fin, it is likely that these morphological (or other 
phenotypic) differences reflect differences in resource use. Further 
experimental studies will provide more evidence. 

5. The sites of sympatry and allopatry do not differ greatly in habitat 
properties and thus food, climate or other environmental features 
affecting the phenotype are likely to be the same. The sites at which 
T. moorii populations live alone are situated at the southern end of 
the distribution area of this Tropheus morph in the same type of 
moderately sloping cobble shore, so it is likely that T polli just 
never successfully colonized this southernmost section. This point 
also includes Grant's (1975) criterion that differences in the zone 
of sympatry should not be predictable from geographic variation 
outside the zone, which clearly is not the case, as it is the same 
color morph of T. moorii. 

6. The sixth criterion requires that independent evidence that similar 
phenotypes actually compete for food should be gained. This 
evidence can be derived from the fact that the > 100 allopatric 
Tropheus populations fill the same function (trophic niche) in the 
complex and species-rich littoral species community. However, 
diet data or stable isotope data are not available so far. 

Taken together, we found two independent lines of evidence for 
adaptive processes underlying body shape differences in closely 
related non-sympatric and sympatric Tropheus populations. Despite 
considerable genetic separation between these populations, all 
sympatric populations occupying deeper water have a relatively 
smaller head, smaller eyes and a more anterior insertion of the 
pectoral fin than do non-sympatric populations. Secondary contact 



of T. polli with T moorii, which have evolved reproductive isolation, 
seems to have enforced a niche shift in both involved taxa, that in 
turn has led to disruptive selection among genetically closely related 
populations and to convergent selection among more distantly 
related populations. The observed concordance in shape in five 
sympatric populations of T moorii, in contrast to the non-sympatric 
populations, together with the shift in resource use, is strong 
evidence of ecological character displacement (Grant, 1972; 
Wassermann and Koepfer, 1977; Simberloff and Boecklen, 1981; 
Schluter, 1986), albeit we are aware that this is not conclusive. Our 
findings are concordant with other cases of early species divergence, 
such as those of Schluter and McPhail (1992), as well as those of 
Adams and Rohlf (2000) on two small and closely related salaman- 
ders (see the review of Losos, 2000), Darwin's finches (Grant and 
Grant, 2006) and nightingales (Reifova et al 2011). 
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